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The search for novel hard materials compared to or 
even harder than diamond, which has the highest mea- 
sured hardness of 96 GPaA and bulk modulus of 443GPa, 2 
has a long history and has stimulated a variety of 
great achievements in high-pressure researchi^^*& Con- 
sequently, many new superhard materials have been pre- 
pared by high-pressure technique, especially after the 
lascr-hcatcd diamond-anvil cells (DACs) was invented. 
In general two groups of materials are powerful candi- 
dates for super-hard materials: (i) strong covalent com- 
pounds formed by light elements, such as polymorphy 
of C3N4, 7 B 6 O f ? and c-BC 2 N. 9 (ii) Partially covalent 
heavy transition metal boride, carbide, nitride and oxide. 
RuOziS and OsB^ii are such examples. Theoretically, 
the nature of hardness has been extensively investigated 
and many new models have been proposed pLiiSiiiiiiiSiiSi 
For the strong covalent materials, hardness can be di- 
rectly derivedfiii^i^*i& while for some metallic transi- 
tion metal-based super-hard materials, it is acknowl- 
edged that bulk modulus or shear modulus can measure 
the hardness in an indirect wayii2*il That is, materials 
with high bulk or shear modulus are likely to be hard 
materials. 

In the present paper, we focused on the bulk modulus, 
mechanical and energetic stability of osmium di-nitride 
(OSN2) with CaF 2 structure^ by calculating the elastic 
constants within the density functional based electronic 
structure method^ For the first time we report that the 
proposed OsN2 compound has very high value of bulk 
modulus (360.7 GPa) which is even higher than that of 
OSO2 with the same structure (347.5 GPa) and is compa- 



rable with that of orthorhombic OsPj2 (365-395 GPaii). 

All first-principles calculations were performed with 
the CASTEP code. 20 The ultrasoft pseudopotential 
(USPP)2i was employed to describe the interaction be- 
tween ions and electrons. Both the local-density approx- 
imation (LDA)2£ and the generalized gradient approxi- 
mation (GGA)2£ were used to describe the exchange and 
correlation potentials. For the Brillouin-zone sampling, 
the Monkhorst-Pack (MP) scheme with a grid of 0.03 A -1 
was adopted i^i The plane- wave cutoff energy is chosen to 
be 550 eV for LDA and 500 eV for GGA calculations. For 
the self-consistent field iterations, the convergence was 
assumed when (i) the total energy difference between the 
last two cycles was less than lxl0~ 6 eV/atom; (ii) the 
maximal force on each atom was below 0.006 eV A" 1 , and 
(iii) the maximal atomic displacement was below 2xl0 -4 
A. We have tested that with even more strict parameters 
the total energy can be converged within 0.002 eV/atom 
for all the systems studied. After getting the equilibrium 
geometry configuration, we applied the so-called "stress- 
strain" method to obtain the elastic constants in that the 
stress can be easily obtained within the density functional 
based electronic structure method* 2 ^ The "stress-strain" 
relation can be described as 

(01, er 2 , 03, 04, er 5 , cr 6 ) = C(ei, e 2 , £3, £4, £5, £e) T ■ (1) 

For the cubic crystal, there are only three non-zero inde- 
pendent symmetry elements (en, C12 and C44). Appling 
two kinds of strains (ei and £4) along the crystallographic 
directions shown in Fig. l(a)-(b), respectively, can give 
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TABLE I: The calculated equilibrium lattice parameters a(A), elastic constants cy (GPa), bulk modulus B (GPa), polycrys- 
talline shear modulus G (GPa), Young's modulus E (GPa), and Poisson's ratio v of typical pure crystals 
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FIG. 1: (Color online) The schematic of strain types: 
(a)ei;(b)£4; (0)53 



stresses relating to these three elastic coefficients, yield- 
ing an efficient method for obtaining elastic constants 
for the cubic system. For the hexagonal crystal, there 
are five independent symmetry elements (en, C12, C13, 
C33 and C44). In order to obtain the additional compo- 
nents C13 and C33, another strain [£3, see Fig. 1(c)] is 
needed. For each strain, in our practical calculations, 
its value is varied from -0.003 to +0.003 with a step 
of 0.0012, then each of three elastic constants takes the 
arithmetic average value of the six steps. The bulk mod- 
ulus is obtained from the elastic constants by the relation 
S=(c u +2c 12 )/3. 

The lattice and elastic constants of diamond, pure plat- 
inum, hexagonal and cubic pure osmium were calculated 
to verify the reliability of the present calculations. It 
is well known that LDA usually underestimates the lat- 
tice constants and overestimates the elastic constants, 
while GGA overestimates the lattice constants and un- 
derestimates the elastic constantsi2Si2i For this reason 
we adopted to use the average of the LDA and GGA re- 



sults as our theoretical estimates. As shown in Table I, 
the theoretical average lattice constant and bulk modu- 
lus of diamond are 3.546 Aand 446.5 GPa, which agree 
well with the experimental values of 3.567 Ai£ and 443 
GPa, 2 with an error of 0.59% and 0.79%, respectively. 
For the cubic osmium and cubic platinum, the theoreti- 
cal bulk modulus values (Pt: 287.2 GPa, Os: 418.8 GPa) 
are also in accordance with other theoretical results (Pt: 
279 GPa^ Os: 417.1 GPa 28 ). With regard to the hexag- 
onal Os, as shown in Table I, the calculated equilibrium 
lattice parameter c/a is 1.581, in good agreement with 
the experimental value of 1.580i22i31 The calculated bulk 
modulus for hexagonal Os is 424.0 GPa, which is also in 
accordance with previous theoretical results (403 GPa, 31 
429.2 GPa 2 ^) and experimental measurements (462±12 
GPa,22 411±6 GPa^ 395±15 GPa 33 ). Based on above- 
mentioned accordance, therefore, we believe that the 
plane-wave ultrasoft pseudopotential (PW-PP) method 
we employed is reliable in investigating the mechanical 
properties of osmium and platinum compounds. 

Now we turn to fully study OsN2 with fluorite struc- 
ture. The results of lattice constant, elastic constants, 
and bulk modulus of OsN2 are listed in Table II. For 
comparison, we have in addition given a calculation on 
platinum dinitridc (PtN2) and osmium dioxide (OSO2), 
and the results are also listed in Table II. Given the 
elastic constants and bulk modulus, the shear modulus 
G and the Young's modulus E can be deduced as fol- 
lows: G =(cn - C12 + 2c 44 )/4, E = 9BG/(3B + G) and 
v=E/ (2G) — 1. These quantities are also shown in Table 
II. 

The key criteria for mechanical stability of a crystal is 
that the strain energy must be positive^ which means 
in a hexagonal crystal that the elastic constants should 
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TABLE II: The calculated equilibrium lattice parameters a(A), elastic constants cy (GPa), bulk modulus B (GPa), polycrys- 
talline shear modulus G (GPa), Young's modulus E (GPa), and Poisson's ratio v of some fluorite and pyrite crystals 
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Os0 2 (fluorite) 
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721.3 
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4.861 
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Ave. 


4.816 


676.8 
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182.8 


347.5 
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579.3 
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87.4 


232.6 
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463.9 
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129 


102 


297.8 (278 d ) 


211.3 


512.7 





.21 




GGA e 


4.862 


668 


99 
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267.3 
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4.819 
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91.8 


288.6 
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103.4 


283.1 





.37 


OsN 2 (pyrite) 


GGA 


4.925 


523 


107 


213 


316 


131 


345.3 





.32 



"Reference 35. "Reference 36. c Reference 26, 37. "Reference 42. e Reference 41. 



satisfy the following inequalities, 

c 44 > 0,cn > |c 12 |, (en + ci 2 )c 33 > 2c 2 3 , (2) 
while for a cubic crystal, 



c 44 > O.cii > |ci 2 |,cii + 2ci2 > 0. 



(3) 



It is straightforward to verify from Table I that the elas- 
tic constants of the hexagonal osmium satisfy formula 
@, implying the stability of hep Os, which is consistent 
with the experimental observation. In the same man- 
ner, from our calculation results in Table I, one can find 
that PtN 2 , OsN 2 and Os0 2 with fluorite structure are 
also mechanically stable since their elastic constants fit 
well in formula Q). The stability of these three crystals 
can also be confirmed by providing the Possion's ratio, 
whose value is usually between -1 and 0.5, corresponding 
to the lower and upper limit where the materials do not 
change their shapes. Note that the present result of bulk 
modulus of PtN 2 is 295.2 GPa. The previous FLAPW 
(full-potential linearized augmented plane waves) calcu- 
lation gives 290 GPa. 26 Remarkably, the two approaches 
agree well, suggesting again the reliability of PW-PP 
method in exploring the structural properties of transi- 
tion metal compounds. On the other hand, we obtained 
the bulk modulus of Os0 2 to be 347.5 GPa, which is 
~13% smaller than that obtained from full-potential lin- 
ear muffin-tin orbital (FP-LMTO) method-S*^ This dif- 
ference may come from different density functional based 
electronic structure method. It reveals in Table I that 
OsN 2 has the highest bulk modulus of 360.7 GPa in our 
series of calculations, this value of OsN 2 is much higher 
than that of other noble metal di-nitride (347 GPa for 
IrN 2 , 190 GPa for AgN 2 , and 222 GPa for AuN 2 37 ). The 
shear modulus of OsN 2 is calculated to be 103.4 GPa, 




FIG. 2: (Color online) Total electron density of OsN 2 at the 
(110) plane. The density at three types symmetry points 
(they are labeled with filled circles, open squares and open 
circles) are approximately 0.8, 0.3 and e A -3 . 



comparable with that of PtN 2 (104.6 GPa) as shown in 
Table I, but much smaller than those of diamond (531.5 
GPa) and Os0 2 (237.1 GPa). Thus compared to dia- 
mond or Os0 2 , OsN 2 cannot withstand the shear stress 
to a large extent. It is interesting to point out that OsN 2 
was implicitly referred in Refill to be unstable or has a 
little bulk modulus. This inconsistency may come from 
the intrinsic technique flaws in WIEN2K code for obtain- 
ing the elastic constants, in which the implanted rhombo- 
hedral distortion is not volume-conservative. When ap- 
plying a non-volume-conservative strain, the expression 
4>eiast = ^2~ c ij £ i £ j f° r elastic energy is no longer accurate, 
resulting in unexpected results in some special cases2&. 

Furthermore, the electronic structure and chemical 
bonding of OsN 2 with fluorite structure are studied by 
calculating its total charge density, Mulliken population, 
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and density of state (DOS). In Fig. 2, we plot the total 
electron density in a (110) plane which cut through both 
the Os and N sites. The bonding behavior of OsN 2 can 
be effectively revealed by analyzing the charge density 
data in real space p(r) at three types crystalline sym- 
metry points, as indicated by filled circles, open squares, 
and open circles in Fig. 2. We found that the charge 
density of these three kinds points are about 0.8, 0.3 and 
e A -3 respectively. Thus the charge density maximum 
lies between Os and N atoms, indicating formation of 
strong covalent bonding between them. Combining the 
fact that each N atom occupies the tetrahedral interstitial 
formed by four Os atoms around it, it is not difficult to 
understand that OsN2 has a low compressibility. Table 

TABLE III: The calculated Atomic and bond Mulliken pop- 
ulation analysis of OsN 2 , Os0 2 and PtN 2 . NM1 and NM2 
denote the first and second non-metal atoms, and M denotes 
the metal atom. 
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Bond (e) 






NM1 


NM2 
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NM1-M 


NM2-M 


NM1-NM2 


OsN 2 


-0.54 


-0.54 


1.09 


1.25 


1.25 


-0.7 


Os0 2 


-0.55 


-0.55 


1.10 


0.98 


0.98 


-0.45 


PtN 2 


-0.45 


-0.45 


0.9 


1.08 


1.08 


-0.34 



III shows bond Mulliken population analysis of OsN2, 
Os02, and PtN2. It indicates that for these three kind 
of materials the bonding is formed between metal atom 
and non-metal atom, while a weak-bonding is formed be- 
tween two non-metal atoms. This is compatible to anal- 
ysis of the electron density of OsN2 in Fig. 2. Table III 
also lists the Mulliken atomic population analysis results, 
from which we can sec the total charge transfer from Os 
to N is 1.10, resulting Os in +1.10 charge state and N 
-0.55 charge state. Therefore, the chemical bonding be- 
tween Os and N has some character of ionicity. It shows 
in Table III that the transferred charge in OsN2 is almost 
the same as that in Os02, and is more than that in PtN2. 
Thus we can make a conclusion that the charge transfer 
effect is more influenced by the metal atom rather than 
the non-metal atom. It is interesting to note that the me- 
chanical properties of OsN2 are also very similar to Os02 
rather than PtN2, as revealed in our above discussions. 

The partial DOS is shown in Fig. 3, no energy gap 
near Fermi level is seen, indicating the metallic nature of 
OsN2. At the Fermi level the total DOS is 1.88 states/eV 
formula units. It reveals that from -18 to -12 eV the 
states are mainly N (2s) states with a small contribution 
from Os (5d and 6p). The states above -9 eV mainly 
come from Os 5c? and N 2p orbitals. 

It was recently demonstrated in both theory and exper- 
iment that the synthesized platinum nitride crystallized 
in pyrite structurei 41 i 42 i 43 The pyrite structure (space 
group number 205), which was also observed in the silica 
recently^ is cubic with 12 atoms per primitive cell. For 
pyrite PtN2, both the four Pt atoms and the midpoints 
of the four nitrogen pairs arrange in fee positions and re- 
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FIG. 3: (Color online) Partial densities of states of OsN 2 . 



suit in NaCl-type arrangement. In addition, each pair of 
nitrogen atoms aligns along one of the (111) directions. 
Besides the lattice constant a, the position (u, u, u) of 
N is the only free structural parameter of pyrite PtN 2 . 
Inspired by these advances, we have also performed a se- 
ries of ab initio total energy calculations to find if OsN2 
favors the intriguing pyrite structure as platinum nitride 
does. Figure 4(a) gives the the energy of OsN2 with the 
internal parameter u varied from 0.23 to 0.40. The plot 
reveals that fluorite OsN2 lies at a local minimum, in- 
dicating its metastable nature. The location of lowest 
total energy is at 0.3614 for GGA calculations, corre- 
sponding to the lattice constant of 4.9246 A of the pyrite 
structure. The bond length of nitrogen pairs in pyrite 
OsN 2 is found to be 1.365 A, even smaller than that of 
pyrite PtN2 (1.51 A). The separation of nitrogen pairs 
in the pyrite OsN2 and PtN2 is nearly the same as that 
of a single-bonded cubic-gauche structure of N observed 
recently^ 4 * 4 ^ which means that the nitrogen pairs proba- 
bly consist some characteristics of covalent bonding. The 
calculated formation energies of pyrite and fluorite OsN2 
at ambient pressure are 0.69 eV and 1.10 eV (per for- 
mula unit), which are smaller than those of pyrite PtN2 
(0.72eV) and fluorite PtN 2 (3.5eV) respectively. Figure 
4(b) plots the enthalpy versus pressure for the fluorite 
and pyrite structures. In the overall range of external 
pressure that we have studied, the enthalpy of pyrite 
OsN2 is always lower than that of fluorite OsN2, imply- 
ing that no first-order phase transition might occur at 
zero-temperature between these two structures. In addi- 
tion, the elastic constants of pyrite OsN2 and PtN2 are 
also calculated and listed in Table II, wherein the results 
for PtN2 show good agreement with those of experimen- 
tal and other available theoretical results. The calcu- 
lated elastic constants of pyrite OsN2 satisfy formula J2J. 
Therefore, it is also a mechanically stable crystal struc- 
ture though its bulk modulus is about 10% lower than 
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FIG. 4: (a) Total energy of OsN2 versus the internal parame- 
ter u; (b) Total energy as a function of the volume per formula 
unit for OsN2 with fluorite and pyrite structure. 



that of fluorite OsN2. 



In conclusion, the OsN2 with fluorite structure is first 
reported to be mechanically stable and have a very high 
bulk modulus of 360.7 GPa by the first-principles calcu- 
lations. The electronic and chemical bonding properties 
have been investigated, indicating that the bonding is a 
mixture of covalent and ionic components. It is found 
that the electronic properties of OsN2 is very similar to 
that of PtN2 with the same structure. As a pyrite-type 
PtN2 and a orthorhombic-type OsN2 have been very re- 
cently synthesized under high pressure and high temper- 
ature conditions j&^ML we expect that the OsN2 as well 
as PtN2 with fluorite structure may be experimentally 
prepared in the future. 
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